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ABSTRACT 

The technique of weak-lensing aperture mass densitometry, so called the (- 
statistic, has recently been popular in actual observations for measurement of 
individual cluster mass. It has however been anticipated that the line-of-sight 
projection by foreground and background matter can adversely affect the cluster 
mass determination with not only substantial error dispersion but also a sizable 
positive systematic bias. Additionally, the finite number of background galaxies 
even at a reasonable observing depth can also introduce Poisson noise to the mass 
estimate. In this paper, we quantitatively investigate the degree of errors sepa- 
rately contributed by the two sources to the mass determination of those galaxy 
clusters with M200 > 10^'^ Mq. We find that the aperture mass of ^-statistic 
turns out to be a mass estimator of much reduced systematic bias, due to the 
cancellation by the positively biased local background mass sheet. However, the 
error dispersion of M200 arising from both projection effect and Poisson noise is 
found to be still sizable (40% — 90%), even for the shear-selected, clean sample 
where multiple clusters located within a suitable projected aperture are removed. 
We also investigate how to remedy this large-error problem in weak lensing mea- 
surements, and propose a plausible alternative mass estimator, M(< 6*1000)5 an 
aperture mass measured within about half the virial radius. The aperture mass 
M{< 6'iooo) is free of bias and has a substantially reduced error dispersion, 39% 
for the worst case of high- 2;, low-mass clusters, that can be smaller than the error 
dispersion of M200 as much as a factor 3. 

Subject headings: cosmology: theory — dark matter — galaxies: clusters: general 
— gravitational lensing 
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1. INTRODUCTION 

The mass function of galaxy clusters has long been recognized as the most convenient 
and important indicator for probing the evolution of structure formation, thereby helping 
determine the cosmological parameters. In addition, the cluster halo mass, when combined 
with the cluster gas mass via the Sunyaev-Zel'dovich effect, provides the opportunity for 
probing the cluster baryon fraction (Umetsu, et al. 2005). Conventional techniques, such as 
measuring the velocity dispersion of gravitationally bound galaxies and the X-ray emission 
profile, have long been employed to measure the cluster mass, assuming chistcr galaxies and 
X-ray emitting plasmas to be dynamically relaxed within the cluster gravitational potential. 

On the other hand, the new technique of mass measurement through weak gravitational 
lensing has been gaining popularity in recent years, with the advantage of not having to 
assume the dynamical equilibrium in the cluster (e.g., Umetsu, Tada, & Futamase 1999; 
Bartelmann & Schneider 2001; Schneider 2005). This methodology was first pioneered by 
Tyson et al. (1990). Various refined techniques for the weak-lcnsing mass determination 
were later proposed by several groups (Fahlman et al. 1994; Kaiser 1995; Bartelmann 1995; 
Seitz & Schneider 1996; Squires & Kaiser 1996; Broadhurst, Takada, Umetsu et al. 2005). 

Among these techniques, the (^—statistic was particularly devised to measure the lens 
mass directly from the tangential component of local gravitational image distortions without 
involving a non-local mass reconstruction (Fahlman et al. 1994; Kaiser 1995). Schneider 
(1996) extended the ^-statistic by generalizing its kernel, which allows one to define an 
optimal measure for the detection of mass concentrations, and this aperture mass technique 
has been apphed to deep optical imaging data to search for clusters (e.g., Erben et al. 2000; 
Umetsu & Futamase 2000; Wittman et al. 2001, 2003; Miyazaki et al. 2002; Dahle et al. 
2003; Hetterscheidt et al. 2005; Schirmer et al. 2006). King et al. (2001) investigated 
the cluster mass measurement infiuenced by interior substructures, and found the measured 
mass as accurate as within 10%. Clowe et al. (2004) studied the effect of asphericity on the 
cluster mass determination, and concluded, under the assumption of an NFW profile, that 
the asphericality effect generally changes the mass estimate by 5% to 10%. 

In addition, several authors have compared the weak-lensing mass with the mass deter- 
mined by the galaxy kinematics and X-ray observations. Reblinsky & Bartelmann (1999) 
concluded that the mass estimates using (^-statistic are significantly more accurate than 
those obtained from the galaxy kinematics. Ettori & Lombardi (2003) studied the mass 
distribution of the rich cluster MS 1008.1-1224 a,t z — 0.302 based on Chandra X-ray and 
FORSl-VLT multicolor-imaging data, and they found that the two mass profiles obtained 
from X-ray and weak-lensing analyses up to 550h~^ kpc are consistent with each other within 
la uncertainty. 
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Irgens et al. (2002), assuming a singular isothermal sphere model for the cluster mass 
profile, compared spectroscopic velocity dispersions, cXp, of 13 X-ray luminous clusters around 
z ~ 0.3 with (TwL of these clusters determined by weak-lensing tangential shear measurements 
out to the cluster virial radius. It was found that among all, two clusters are in strong 
discrepancy, with (Jwl > 2o"p and (Jwl ~ 2ap, whereas the rest are in fair agreement, with 
(cp/cwl) ~ 1- Though these exceptional clusters may be in dynamical non-equilibrium, 
another possibility may arise from the projection effects of other mass concentrations (Cen 
1997; Rebhnsky & Bartelmann 1999; White, van Waerbeke, & mackey 2002; Padmanabhan, 
Seljak, & Pen 2003; Hamana, Takada, & Yoshida 2004; Henawi & Spergel 2005) and/or local 
filamentary structures (Metzler et al. 1999; Metzler, White, & Loken 2001) along the line-of- 
sight. Using A'^-body simulations White et al. (2002) studied the completeness and efficiency 
of weak-lensing cluster surveys on the basis of their mass-selected mock cluster sample and 
found that the line-of-sight projection effects can be quite serious due to the broad lensing 
kernel. In the cluster mass estimate based on the convergence map, they found a positive 
bias of ~ 20 — 30% with a substantially larger error dispersion that can even occasionally 
yield negative lens masses. Metzler et al. (1999; 2001) studied the projection effects on 
weak-lensing mass estimates for massive clusters caused by the local large-scale filamentary 
structures. Including the projection effects from local matter within a sphere of 128/i~^Mpc 
radius, they found the lensing convergence maps to yield an positive mass bias of 30% 
and the mass error dispersion of ~ 0.3 for massive clusters at a redshift of z — 0.5. In fact, 
these problems of weak-lensing mass determination, i.e., positive mass bias and large mass 
error dispersion, have been alluded in earlier works (Cen 1997; Reblinsky & Bartelmann 
1999). Further, the cluster halo triaxiality itself can cause a bias in the lensing-based mass 
estimation (Clowe et al. 2004; Hamana et al. 2004; Oguri et al. 2005), while it is hkely to 
have less effect on the X-ray cluster mass estimate (Gavazzi 2005). Thus, despite that weak 
lensing offers a unique tool for the measurement of cluster masses without any assumption 
of their equilibrium state, it can however suffer from the projection effects. Such problems 
are less significant in X-ray or spectroscopic velocity-dispersion measurements. 

The present study aims to investigate the errors in weak lensing cluster mass measure- 
ments as well. However, this work differs from the aforementioned previous works, in that 
we attempt to simulate the actual wide-field weak lensing measurements, with numerically 
simulated shear data as closely resembling the observing data as possible. In particular, we 
shall focus on the bias errors and random errors pertinent to local weak-lensing measure- 
ment of (^-statistic. Moreover, at a given observing depth, the finite number of background 
galaxies can introduce non-negligible Poisson noise convolved with the projection error in 
the measured data. We shall quantify, in this work, the regime for which the projection 
effect dominates, and the other regime where the Poisson noise dominates. 
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This paper is organized as follows. In Sect. 2 we describe our cosmological A'"-body 
simulations, weak lensing simulations, and the construction of a mock cluster catalog. Details 
of our shear-based mass estimator and mock observations are presented in Sect. 3. In Sect. 4 
we apply the shear-based mass estimator to our simulated weak lensing observations, and 
examine the statistical properties of the errors in weak lensing cluster mass estimates. The 
radial mass error profiles are discussed in Sect. 5. Based on the radial error profiles, we 
propose a plausible alternative cluster mass indicator that has much reduced error dispersion. 
We finally present the discussions and a summary in Sect. 6. 

2. Simulations 

2.1. Particle A^-Body Simulation and Weak Lensing Simulation 

We use the Gadget code (Springel et al. 2001) to run 10 sets of independent simulations 
with 128^ dark matter particles for a ACDM model in a lOOh"^ Mpc cubic box. The cosmo- 
logical parameters are flm = 0.3, fl\ = 0.7, the Hubble parameter h = i^o/lOOkm/s/Mpc = 
0.7, and the linear mass fluctuation amplitude erg = 0.94. The mass resolution is 3.2 x IO^^Mq 
per particle. We also conduct one simulation of 8 times higher mass resolution with the same 
initial condition for checking the convergence of the result. 

General methods of backward ray-tracing have been detailed in Wambsganss et al. 
(1998) and Jain, Seljak, & White (2000). In the present study, on the other hand, we adopt 
a simpler linear approximation for the investigation of weak gravitational lensing. The shear 
matrices are generated on every lens plane with 4096^ resolution, and we trace uniform 1024^ 
rays backwards from the observer plane {z = 0) to the source plane {z = zs)- In the weak 
lensing approximation that employs the Born approximation, each photon travels along the 
un-perturbed trajectory, and the accumulated shear is the linear combination of gravitational 
shear given by every lens plane. Throughout this paper, we adopt a single source plane 
located at ^5 = 1, which is the typical value for the mean redshift of background galaxies 
in an actual weak lensing survey with limiting magnitude R ~ 25.5 mag (see Hamana 
et al. 2004). We then randomly choose a Ideg^ field of view from the large projected 
simulation cubic box, and generate 36 lensing maps of Ideg^ from each simulation. Ten sets 
of independent lOOMpc cosmological simulations were used to avoid the same clusters to be 
observed repetitively. Hence, the effective survey area is 36 X 10 X Ideg^ = 360deg^ 

To make sure the difference between 128'^-particle simulation and 256'^-particle simu- 
lation indeed negligible for cluster mass measurement, we compare 256'^-particle simulation 
with 128'^-particle simulation of the same initial condition for the weak lensing estimation of 
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cluster mass. It is found that two results are almost indistinguishable, for the following rea- 
son: The resolution of cluster mass estimation is limited by the finite galaxy number density, 
Ug ~ 30 galaxies per arcmin^ which corresponds to the grid resolution of about 300 x 300 
per deg^. On the other hand, even a small cluster of 10^^ Mq of redshift z = 0.2 contains 
3000 particles within a projected area of 10 arcmin^. The ten-times higher dark-matter 
particle number density than the background galaxy number density already warrants little 
dumpiness in mass distribution to arise from the discrete particle effect, thereby ensuring 
correct mass estimation even with lower resolution simulations. 



2.2. Mock Cluster Catalog 

Wc use the fricnds-of-fricnds algorithm (Huchra & Gcllcr 1982) with linking-length 0.2 
to search for clusters. After a cluster is identified, the center of mass position Vc is then 
computed. We define the M200 of a cluster by increasing the enclosing spherical radius 
around the cluster center rc till it satisfies the criterion of r2oo, within which the mean 
interior density is 200 times the mean density p{z) — ^rn{z)pc{z) of the universe at the 
cluster redshift. Similarly, we define the cluster mass enclosed within ta^ by 

Att 

Ma. = -jA^p{z)ri^ (1) 

with Ac being the mean overdensity within ta. with respect to p{z). 

In the end, we project the 3-D r^. to obtain the projected cluster center, 0c, and the 
projected r2oo to 6200, respectively. The peak position in the 2-D lensing map dose not neces- 
sarily coincide the center of mass position determined in 3-D particle simulation, especially 
in the case where the halo in 2-D projection contains substructures or mergers. Prom the 
observational viewpoint, the cluster center can be chosen from the local surface density peak. 
Nevertheless, we expect the difference not to be so serious as to affect the statistics of our 
mass estimate, since the difference between two centers is small compared with the inner 
most radius where the cluster mass estimate is made. 

Through the above procedure, we obtain a mass-selected sample of about 4000 clusters 
with M200 > IO-'^^Mq between lens redshifts z — 0.2 to 0.8 in our effective survey area 
360deg^. Repeating clusters do occur in our randomly rendered maps but only at different 
cluster orientations and with different background and foreground lenses. For the purpose 
of the present work, which focuses on the measurement errors, repeating clusters can be 
regarded as independent samples since the noise is different. 
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3. 



Method 



Our method for cluster mass estimation is based on the gravitational shear field, which 
is different from the method based on the convergence field used by Metzler et al. (2001) 
and White et al. (2002). Using the tangential component of galaxy shear for a lensing mass 
determination is often a preferred one, since it makes direct use of the locally observed galaxy 
ellipticities around the target cluster and is widely adopted in recent observations (Clowe 
et al. 2006; Bardeau et al. 2005; Jee et al. 2005). By contrast, the convergence-based 
method involves non-local mass reconstruction from the observation of galaxy ellipticities 
over a large field. Mathematically, the two methods are equivalent, and both exhibit the 
mass-sheet degeneracy. But in practice the local shear-field measurement has one additional 
degree of freedom in choosing an appropriate nearby control field as the background mass 
sheet to be subtracted. 



We begin by briefly summarizing the basic formalism of weak gravitational lensing. A 
general review of weak lensing can be found in Bartelmann & Schneider (2001). 

The deformation of an infinitesimal ray bundle due to gravitational deflection is de- 
scribed by a mapping between the two-dimensional source plane and the image plane as 



where /3 and denote the angular position on the source and the image plane, respectively, 
and Aij is the 2x2 Jacobian matrix of the lens equation. In the weak lensing regime, the 
Jacobian matrix Aij is symmetric and can be decomposed as 



where is the shear matrix deflned as — 71(73 + 72(71 with 7^ being the components of 
complex gravitational shear 7 := 71 -|- 272, (7j being the 2x2 Pauli matrices, and k, being 
the lensing convergence responsible for the trace-part of the Jacobian matrix. The lensing 
convergence k is a line-of-sight projection of the matter density contrast 6 = (p^ — p)/p out 
to the source plane {S) weighted by certain combination g of comoving angular-diameter 
distances (e.g., Jain et al. 2000), 



3.1. Weak Lensing Formalism 



(2) 




(3) 
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where a is the cosmic scale factor, and x is the co-moving distance; is the matter surface 
density Ej„ = Jdx ci{Pm — p) with respect to the cosmic mean density, and E^it is the critical 
surface mass density of gravitational lensing, 

Ee. = (5) 

with Di, D5, and Di,s being the angular-diameter distances from the observer to the lens, 
from the observer to the source, and from the lens to the source, respectively. The gravi- 
tational shear field is related with the lensing convergence field in a non-local manner (e.g., 
Bartelmann & Schneider 2001). The relation between k and 7 is expressed in Fourier space 
as (Kaiser & Squires 1993) 



Ml) = f^^^ii)^ (6) 

2/1/2 
II + 11 

(8) 



72(0 = TF-^'^IO, (7) 



where k{l) is the Fourier transform of k,{0), is the Fourier transform of ^i{0), and I is 
the Fourier variable conjugate to the angular position on the sky. Equations (6) and (7) 
can be used to invert the gravitational shear field to the lensing convergence field. 



3.2. The C- Statistic Mass Estimator 

Jee et al. (2005) performed a weak lensing analysis on a 2; = 0.8 cluster, and found the 
mass estimate obtained from the aperture densitometry (Fahlman et al. 1994), or so-called 
the (^-statistic, to be very close to that obtained from the nonlinear iteration-convolution 
method of Kaiser & Squires (1993). They also investigated the (^—statistic within a projected 
aperture of 0.15r2oo, where the weak lensing assumption almost breaks down, and showed 
only about 10% mass error. Although the ^-statistic can be accurate from the linear regime 
to the weakly nonlinear regime, this method however cannot avoid the projection effect 
arising from the foreground and background matter along the line-of-sight. In what follows, 
the mass errors caused by projection as well as other error sources will be quantified. 

The observed image distortion of background galaxies can be directly used to derive 
the projected gravitational mass of clusters. The aperture mass estimate within the angular 
radius 9i, M^{< 9i), in terms of the tangential component 7^ of gravitational shear can be 
expressed as 

M^(< ^1) = 7r(^iL'i)2E,HtC(^i), (9) 
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using the C-statistic defined as 

C(^i) -.^ ^_l,j^, jyin9{M0))^n< 9i)-m,02), (10) 

where (...) denotes the azimuthal average, and K, = Sjn/Scrit is the mean convergence. Equa- 
tions (9) and (10) show that the cluster mass can be measured from the galaxies ellipticity 
within the annulus bounded by 9i and 02 located just outside the mass to be measured. 

As revealed in Eq. (10), the ^—statistic yields the mean convergence interior to 6i, 
subtracted by the mean background within the annulus between 9i and 02, k(< 9i)—K{9i, 6*2). 
Hence, as long as k{9i,92) <^ 'k{< 9i), the enclosed mass within 9i can be obtained by 
multiplying ( by the area 7rEcrit(-Di:,^i)^. As the inner radius 9i can almost be arbitrarily 
chosen to obtain the aperture mass within 9i, when 9i is chosen to be ^200, the aperture mass 

is approximately the cluster mass, M2oo- Obviously, the aperture mass is smaller than 
the enclosed mass by a negative compensating mass that serves to remove the contribution 
from a background uniform mass sheet, and the degree of deviation depends on how steep the 
density profile is. In §4.3, we will demonstrate empirically that the negative compensating 
mass can actually correct for the systematic positive bias resulting from the projection effect. 

The (^-statistic is a spherically symmetric mass estimator. It gives rise to some errors 
for the typically irregular cluster. Nevertheless, this effect has been estimated less than 10% 
(Clowe et al. 2004). To avoid this effect to inffuence our cluster mass estimation, the number 
of clusters in the sample is crucial. Our sample of 4000 clusters contains clusters of different 
sizes and shapes, and the effect of non-spherical lens is expected to average out. Moreover, 
to avoid strong lensing, our mass estimator avoids the shear measurement near the cluster 
center (see §4.4). 

The choice of the parameter 92 may also affect the cluster mass estimate. For example, 
a small ^2 will generate large Poisson noise since the galaxy number for shear estimate within 
the annulus bound by 9i and ^2 is not sufficiently large. On the other hand, if ^2 is too large, 
the cluster mass measurement can be contaminated by the neighboring clusters in a crowded 
field. In practice, we nevertheless obtain similar results for a = ^2/^1 = 1.15,1.2,1.3 and 
1.4 with 9i = 9200- Bardeau et al. (2005) and Jee et al. (2005) adopted the parameter a of 
1.22 and 1.15, respectively. Throughout this paper, unless otherwise noted, we use a — 1.2 
for the measurements of M200. 

For a projected lens system, (7^) is produced not only by the cluster itself but also by 
the projected neighboring clusters and large scale structures. As extreme examples. Figs. 1 
and 2 show two convergence maps that give rise to large errors in the cluster mass estimate. 

Figure 1 shows a simulated cluster of underestimated M200, since a high mass cluster of 
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different redshift is located just in the annulus bounded by 9i and ^2- The underestimated 

mass can become negative. Figure 2 shows a cluster of overestimated M200, as another cluster 
of different redshift is enclosed by 61. In this case the cluster mass is overestimated by a 
factor >2. Both Figs. 1 and 2 have fields of view of 1. 2^200 defined by the central objects. 
Peculiar clusters of this kind are included in our raw sample of 4000 clusters. 

3.3. Adding Noise 

Another important factor that strongly influences the mass measurement is the galaxy 
intrinsic ellipticity, which plays a role as random noise in the shear measurement. In the 
weak lensing regime, the shear induced by gravitational Icnsing is some slight distortion of 
galaxy images, e.g., 10% distortion for cluster lenses. However, individual galaxies have an 
intrinsic ellipticity of random orientation with considerable dispersion, about 30 — 40%. Such 
orientation noise can only be reduced by taking average over some sufficiently large projected 
area around the lens. 

In our simulation, the two components of the complex galaxy intrinsic ellipticity, ej^* 
{a = 1,2), are generated by Gaussian random numbers with dispersion a^/\/2 per compo- 
nent. In the weak lensing approximation where k <^ 1 and I7I -C 1, the net ellipticity e"*^* is 
a linear combination of galaxy intrinsic ellipticity e^^^ and gravitational shear 7^, that is, 

C = C* + 7a. (11) 

In what follows, we assume a weak lensing survey of mean galaxy surface number density, 
Ug — 30 arcmin"^, and intrinsic ellipticity dispersion, cr-y = 0.4, in our simulations. 

3.4. Mass and Redshift Bins of Cluster Sample 

Given a finite number of sample clusters in the simulation, we focus on those clusters 
that can produce sufficiently strong signals for the weak-lensing mass determination. This 
issue is important in the context of this study, as we are interested in the observational 
uncertainties in the cluster mass determination. Clusters of unsuitable mass and redshift 
ranges will yield unreasonably high uncertainties to provide useful astrophysics information. 
Our sample of clusters is selected based on the mass (e.g.. White et al. 2002), whereas 
practical weak lensing samples of clusters are "shear- selected" , that is, weighted by both the 
mass and the redshift (e.g., Reblinsky & Bartelmann 1999; Hamana ct al. 2004). A concise 
overview of how clusters of viable mass redshift ranges are determined is given below. It 
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provides a criterion for us to select those "clean" clusters from 4000 clusters in the raw- 
sample. 

Typical cluster search schemes use the convergence map instead of the shear map (e.g., 
White et al. 2002; Padmanabhan et al. 2003; Hamana et al. 2004). In this method one first 
adopts a Gaussian filter to smooth the raw convergence k map, followed by identification of 
local peaks in the smoothed convergence map with a certain threshold. The threshold signal- 
to-noise ratio (S/N) in the smoothed K-map is given as u = Kg/uk with rms noise in the 

KG-map, (Jk = o"7/(4vr^^^n3)^/^ with the Gaussian smoothing scale Oq = FWHM/ y/ 4 ln(2) ~ 
0.02((T^/0.4)(np/30arcmin-2)-V2(^g/l')-i. 

An optimal choice for the threshold v depends on the concentration parameter of the 
NFW profile, Cnfw! the cluster mass M^i and the lens redshift z^. Moreover, the optimal 
threshold v can be quantified in terms of the survey completeness and efficiency. For a weak 
lensing survey with Ug = 30 arcmin"^ and zs = I, corresponding to an actual survey with 
limiting magnitude R ~ 25.5 mag under a subarcsec seeing condition, Hamana et al. (2004) 
found that the threshold i/ ~ 4 with ~ 1' provides an optimal balance between survey 
completeness and efficiency. At the optimal threshold i/ ~ 4, the lower mass detection limits 
are then = IO^^Mq, 2 x lO^^M©, and 4 x lO^^M© for lens redshift ranges of 0.2 < zl < 0.4, 
0.2 < Zl < 0.6, and 0.2 < zl < 0.8, respectively. We therefore divide the cluster mass into 
three mass bins: 

IO^^Mq < Ml < 2 X IO^^Mq, (12) 
2 X IO^^Mq < M2 < 4 X lO^^Mo, (13) 
M3 > 4 X IO^^Mq, (14) 

where Mi, M2, and M3 are defined with the M2oo- The above three redshift bins are also 
adopted in this study. The cluster numbers in different mass bins and different redshift 
intervals in our sample are given in Table 1. 



4. Mass Probability Distribution Function 

For each cluster in our mock cluster sample, we determine the cluster mass AftruQ using 
the A^-body simulation data and the lensing mass Miens = using the ^-statistic (10). Wc 
note that both M^rue and Miens are projected masses enclosed by a cylinder of certain radius 
Dl9 which is relevant to gravitational lensing. For the true mass Mtme, the integral along 
the line-of-sight is taken only within r2oo from the cluster center, while it is taken over the 
hne-of-sight from 2; = to 2; = ^5. When the angular radius 9 is taken to be ^200 = ^200 1 Dl, 
then M(< ^200) ~ -^200- Collecting all clusters of M200 > IO^^Mq from the simulation data. 
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we then construct the observed mass probabihty function with 

However, as will be elucidated later, a substantial fraction of lenses are suffered from 
the projection effect as well as strong observational noise due to intrinsic galaxy ellipticities, 
both of which yield considerable errors in the mass measurement. To reduce such errors 
in the mass determination, we conduct a simple-minded "clean sample" procedure, which 
removes those clusters from our raw sample that obviously contain nearby objects in the 
convergence map, and/or those clusters whose mass is below a redshift-dependent limiting 
mass (§4.2). 

We nevertheless note that even in the clean sample, mass estimates from the C-statistic 
still contain residual mass errors due to local filamentary structures around the target cluster 
as well as intervening large-scale fluctuations (i.e., cosmic shear) projected along the line-of- 
sight, as shown in Fig. 3. The details of the clean procedure will be presented in Sec. 4.2. 
The numbers of clean clusters in three mass bins and three redshift intervals are listed in 
Table 2. Most of the removed samples contain contaminant objects of different redshifts due 
to chance projections along the line-of- sight. Only a small fraction contain physically nearby 
objects and mergers, and the parenthesis in Table 2 gives the percentage of removed samples 
that belonging to this category. 



4.1. Raw Sample 

The observational errors of M200 is examined from our raw sample. The result derived 
from the raw sample serves as a baseline for comparison with those derived from the clean 
sample. The mass probability distribution function is defined as -P(/i) = dn/ {Nzdii), where 
Nz is the total cluster number in a given redshift bin, d/i the mass error bin, and dn the 
cluster number in the mass error bin at a given redshift. 

Two cases are separately investigated: the noise-free case that yields Fig. 4, and the 
noisy case that yields Fig. 5. The noise- free case contains measurement errors due solely 
to the projection effect, and the mass probability function is positively skewed with a tail. 
Nevertheless, we find the peak at // 1, which is in agreement with White et al. (2002). 

On the other hand, the noisy case contains contaminations from both projection and 
galaxy ellipticity effects. The random Gaussian noise due to intrinsic galaxy ellipticities 
adds substantial errors to the gravitational shear estimate, giving rise to a mass probability 
function considerably more symmetric and broader than the noise-free case. 
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The error bars in both Figs. 4 and 5 represent the intervals of 68% confidence. We note 
that for the noisy case, the confidence intervals are obtained by averaging over 20 independent 
galaxy realizations. Due to the Gaussian nature of the galaxy ellipticity noise, the S/N can 
be improved in a predictable manner as ^/N^, where A*"^ is the number of background galaxies 
used to determine the cluster mass. Hence, high-redshift clusters, being smaller in angular 
size, have noisier mass estimates than low redshift clusters. Furthermore, as high-redshift 
clusters are not optimally located for the zs — 1 source galaxies, only massive ones can yield 
reasonably acceptable mass determination. 

The mass-error dispersion of the respective mass probability functions in three different 
redshift intervals is given in Table 3. It is clear that the error dispersion for the raw sample 
is so large that the mass measurement can hardly be in good use for construction of the mass 
function. Again, Table 3 serves as the baseline for comparison against the error dispersion 
resulting from the clean sample. 

4.2. Clean Sample 

We apply a "clean" procedure to the raw cluster sample (§4.1) in order to remove clusters 
that would yield unreliable mass estimates from the ^-statistic mass estimator (§3.1). Our 
clean procedure consists of the following two steps: The first step is to remove clusters whose 
masses are below a weak-lensing detection limit. Adopting the result of Hamana et al. (2004) , 
we first set lower limiting masses of M200 = 4 x 10^^ M©, 2 x lO^^M©, and IO^'^Mq, for clusters 
in redshift ranges of 0.2 < 2; < 0.8, 0.2 < 2; < 0.6, and 0.2 < ^ < 0.4, respectively (see §3.3). 
The clean sample contains only these clusters selected by the redshift-dependent limiting 
masses M2oo(-2l) above. The sample after this first cleaning procedure can be regarded as 
an effective "shear-selected" sample of clusters. The second step is to remove those clusters 
which have detectable "nearby" clusters in the projected map located within 1. 2^200 from the 
target cluster. After the above removal procedure the remaining randomly rendered maps 
contain only the clean clusters that constitute our clean sample. Table 2 shows the cluster 
numbers of the clean sample in different mass and redshift bins, as well as the percentage of 
removed clusters that have genuinely nearby contaminants in the physical space. 

The clean procedure reduces the projection effect in the shear-based cluster mass es- 
timate, and the benefit of "clean" is found to be only moderate as revealed by the error 
dispersions of the clean sample in Table 4. A comparison between Tables 3 and 4 shows that 
the clean sample has 15% improvement for the least massive clusters of low-redshift to to 
50% for the most massive clusters of high-redshift in the noise-free error dispersion. When 
observational noise due to intrinsic galaxy ellipticities is included, the clean sample is at 
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most 30% better in error dispersion. Such a result cannot be considered to be satisfactory. 
Clearly, there is still much room for improvement in error reduction. This will be the subject 
of §4.4 and 4.5. 

4.3. Moments of Clean-Sample Mass Probability Distribution P\mctions 

Figure 6 depicts the mass probability distribution functions P(/X2oo) of //200 = -^true(< 
^2oo)/-^iens(< ^200 ) in uoisy and noise-free simulations for different mass-bins and redshift- 
intervals derived from the clean sample. The statistical properties of a probability function 
can be characterized by its distribution moments, such as the mean, the dispersion, the skew- 
ness, and other higher-order moments. We compute these first three moments to quantity 
the mass probability function for the 3 mass bins and 3 redshift intervals. The results are 
listed in Tables 4, 5 and 6. 

The mass probability function tends to have a positive tail and be asymmetric. The pos- 
itive offset can be characterized by the skewness S3 tabulated in Table 6. The skewness is de- 
fined to be the third moment of the distribution function: S3 := '^f^iifJ'i — fiY / {{N — l)o-^)- 
This parameter is a measure of non-Gaussianity originated from the projection effect, since 
the contaminant lenses are non-Gaussian. However, inclusion of the galaxy ellipticity noise 
substantially reduces the skewness, thus making the mass probability function more sym- 
metrical about the peak. 

Note from Table 5 that the mean is found to be only slightly greater than the peak 1) 
shown in Fig. 6 as a result of the small positive skewness of the distribution. To be specific, 
the mean has a positive bias of less than 10% for clusters of all mass and redshift ranges 
in our clean sample. Moreover, the positive bias increases only slightly with increasing lens 
redshift. The fact that the ^-statistic mass estimator has a small bias is surprising, since it 
is at variance with a previous study. White et al. (2002) adopted to use the convergence 
map for a weak-lensing mass measurement, and showed a 20% positive mass bias for massive 
clusters of intermediate redshifts. Such a positive bias results from the projection effect 
by local large-scale structures surrounding the target galaxy cluster that accounts for an 
additional 20%-30% lensing strength (Metzler et al. 1999; Metzler et al. 2001). 

In contrast to the /t-based mass measurement, the local shear-based mass measurement 
with ^-statistic turns out to have a compensating effect to neutralize the anticipated positive 
bias. Weak lensing measurements are subject to the mass-sheet degeneracy (Bradac et al. 
2004), which renders the uniform background matter undetectable. In the K-based mass 
measurement, the overall field tends to be so large as to permit an accurate global inversion 
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of the convergence field. The background so defined is therefore the global background over 
a large field. On the other hand, the shear-based measurement is a local measurement and 
the background matter density is defined only locally in the projection space. The second 
term in Eq. (10) is exactly the local background surface density defined over a control field 
within an annulus just outside the target field. Hence while the environment surrounding the 
cluster is rich in matter, the matter density in the control field just outside the aperture of the 
target field can be correspondingly high as to cancel the enhanced background contribution 
along the hne-of-sight of the target field. The cancellation significantly reduces the positive 
bias of the lensing mass. 

To show the cancellation of the bias, we compute the ratio of the noise-frcc local mass 
AM(6'i,6'2) = 7rScrit-D|(6'| — 61)k{6i,62), contained in the second term in Eq. (10), to the 
true cluster mass Mtrue, where all matter outside 50/i~^Mpc from the target cluster has been 
excluded for examination of the local contribution to the projection effect. Here we take 
01 — 0200 and 02 = l.29i. Table 7 shows the means and the variances of this ratio for all 
clusters in the clean sample. The mean ratio ranges from 18% to 26% with a considerably 
small variance. 

Indeed, the mean ratio has a correct value to largely cancel the well-known mass over- 
estimation of 20 — 30% produced by the local projection effect (Mctzlcr et al. 1999; Metzler 
et al. 2001; White ct al. 2002). This result quantitatively demonstrates the self-cancellation 
of two environmental contributions to be at work with (^-statistic. The possibility of such a 
cancellation was previously alluded by Mctzlcr ct al (1999), but the small number of sample 
clusters prevented them from drawing a quantitative conclusion. 

Coming back to Table 4, we find the error dispersion for the noise-free case tends to 
be small at redshift intervals in between 0.4 < Zl < 0.6. This is because the source galaxy 
redshift Zs = 1 favors lensing signals from this lens redshift. The noise-free result for the 
high mass bin, M3, at 0.4 < zl < 0.6 can be compared with the result of Metzler et al. 
(2001) for z — 0.5. They obtained an error dispersion a — 0.26 in the ratio between the 
estimated M200 and the true M200, which is in good agreement with the noise-free error 
dispersion in our clean sample, a — 0.28. On the other hand, for the noisy case, the smallest 
error dispersion occurs at lens redshift interval 0.2 < zl < 0.4. It means that the errors 
contributed by galaxy ellipticity noise is significant, and the larger angular size of low-z lens, 
having a larger sampling annulus, helps reduce the sampling noise. Table 4 also shows that 
the cluster mass measurement errors can be significantly enhanced by the galaxy ellipticity 
noise by more than almost 50% for all lenses, except for the highest mass and lowest redshift 
bin. This indicates that the dominant error source in the weak lensing M200 can be the galaxy 
ellipticity. This finding motivates us to proceed on seeking an alternative mass estimator 



- 15 - 



below. 



4.4. Radial Mass Error Profile 

We turn to the radial profile iJ,{6) = Miens (< ^)/Mtrue(< measured for individual 
clusters of the clean sample. We fix the outer radius 62 to be 1. 2^200 and vary the mass- 
aperture radius 9i. The results are presented in Fig. 7 for the noise-free and noisy cases. 
The innermost radius is chosen to be 0. 2^200 so as to stay in the observationally weak-lensing 
regime. For a typical cluster in our clean sample, r2oo ~ l/i~^Mpc and the innermost radius 
limit is about 200/i~^kpc, which is mostly outside the scaling radius of the NFW profile. 

The error bars in Fig. 7 represent 68% confidence levels around the estimated mean 
Ijl{9), obtained over 20 independent background reahzations similar to Fig. 5. The deviation 
of the mean-/i from unity indicates a systematic bias. Very little systematic bias is detected 
interior to ^200 in Fig- 7. 

The error bars in Fig. 7 are seen to decrease with decreasing radius. In particular, the 
mass error caused by intrinsic galaxy ellipticitics (i.e., the difference between error bars of 
noisy and noise-free cases) decreases mostly noticeably for ^200 > > 0. 7^200 • This feature 
arises from the peculiar feature of (^-statistic, where the enclosed mass within a smaller radius 
is measured by a larger annulus. Such an error stays relatively constant for 9 ^ 0.7^200i which 
can be understood from the geometry. The annulus area enhancement begins to saturate 
when the inner radius becomes smaller than half the outer radius. On the other hand, 
the projection errors continue to decrease at the ever decreasing radius. This trend can be 
understood from the second equahty of Eq. (10). For a small 9i, the interior convergence 
k{9i) is much greater than the average convergence in the annulus k{9i, ^2), which contains 
lensing contributions of other redshifts. It therefore becomes increasingly insensitive to most 
contaminant lenses for a decreasing 9i . 

4.5. Alternative Weak-Lensing Mass Estimators? 

Section 4.3 already reveals the considerably large errors arising from galaxy ellipticity 
in the weak lensing mass measurements. On the other hand. Sec. 4.4 also shows that despite 
a sizable error in the measured M200, the error becomes substantially small for the fractional 
mass well within r2oo- This tendency has also been noticed by Metzler et al. (2001) and Cen 
(1997) previously for noise-free mass estimations. 

In order for the weak-lensing measurement to yield a reliable cluster mass estimate. 
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we seek an alternative mass estimator that is contaminated by the eUipticity noise and the 
projection effect to a much lesser degree. In Table 8, we enlist the ever-decreasing mass 
errors for the 2D projected mass measured over the decreasing radii, ^500, ^looo; and ^^1500- 
Here 9/^^ = r^^/Di is the projected angular radius of the 3D interior mass Ma^ defined by 
Eq. (1). These radii correspond, on average, to = 0.76'2oo, 0.5^^2oo, and 0.4^^200) respectively, 
in weak dependence on the lens mass and redshift. Projected masses evaluated at these radii 
contain significantly smaller errors than M200, and may serve as alternative mass indicators. 

A similar proposal was put forth for the mass estimator of X-ray clusters. With cos- 
mological simulations, Evrard ct al. (1996) found that X-ray mass estimates arc remarkably 
accurate when evaluated at radii from rsoo to r25oo- Recent X-ray observations measure the 
cluster mass at rsoo (Finoguenov et al. 2002) and at r25oo (Allen et al. 2002). Moreover, Jee 
et al. (2005a) compared the cluster mass profiles derived from X-ray and from weak lensing, 
and found increasing agreement with decreasing radii. The radius r25oo is about 0.3r2oo for 
our sample clusters. In fact, the mass within 0. 8^200 may be accurately measured by weak 
lensing even at a lens redshift of zl ~ 0.8, as discussed by Jee et al. (2005b). However, the 
cluster mass within 0.3r2oo can be sensitive to the complex interior structures of a cluster, for 
example, containing two cores in a cluster recently undergoing a merger, and/or the baryonic 
physical processes such as radiative cooling and galaxy formation. Moreover, for a dual-core 
cluster, r25oo can be ambiguous to define. 

We therefore suggest that M(< ^1000) be a better weak-lensing mass estimator. On one 
hand, it has an acceptable worst mass errors (< 39%) for the 3 mass bins and different redshift 
intervals, as revealed in Table 8. On the other hand, the radius is still sufficiently large, 
~ 0.56*200, so that the mass estimator can be insensitive to the cluster interior structures. 
Finally, M(< ^1000) is a significant fraction of M200 and it should still follow the well-known 
similarity scahng of the CDM gravitational collapse. 

5. CONCLUSIONS 

This paper reports a systematic study on the accuracy of cluster mass measurements 
through weak lensing observations. Specifically this work takes into account the mass errors 
introduced by the projection effect and by the Poisson noise of finite number of randomly 
oriented background galaxies. Among these sources of mass errors, the projection effect 
has been reported to yield a non-negligible systematic positive bias (~ 20%) in the cluster 
mass estimate based on weak lensing. Adopting the local shear-field measurement using 
^-statistic, which is often used in actual observations, we nevertheless found that such a 
positive bias can be largely canceled by the positively-biased local background mass sheet. 
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That is, the C-statistic measurement can provide a bias-free cluster mass estimate. 

In this paper, we also report that the error in M200 determination is expected to exceed 
50% even for a moderately deep observation (R ~ 25.5 mag), regardless of lens mass and 
redshift. Even after a clean procedure that removes clusters with detectable companions 
in the projected map, the mass error is still substantial, exceeding 40% for a ~ 25.5 
observation. Mass errors can scatter observed data of one mass bin into another bin when 
constructing the cluster mass function, and smear out the mass function. In the mass 
range where the mass function has a large gradient, i.e., M > M^, such mass errors can 
greatly distort the mass function. The measurement error will eventually propagate into the 
determination of cosmological parameters, such as the matter density parameter and the 
matter fluctuation amplitude as, which rely critically on an accurate mass function. The 
so-called self-cahbration was devised to correct for the systematic errors of this kind (Hu 
2003), but the random errors are un- removable even with the self-calibration. 

To significantly reduce the mass error in weak lensing measurements, we also suggest 
the possibility of an alternative lens mass, which has a considerably smaller error. It is 
an interior mass well inside ^200- For example. Fig. 7 shows that the error in M(< ^500) 
(mass interior to approximately 0. 7^200) is already noticeably smaller than M200, the error 
dispersion of M{< ^1000) is at most 39% for all detectable clusters and the error dispersion of 
M{< 61500) at most 32%. Comparing Table 8 with Table 4, we find Miooo to have more than 
a factor 2 in error reduction for all detectable clusters; for low- 2; lenses, the error reduction 
can be as large as a factor 3. Given these results, we suggest that M(< ^1000) be a better 
mass variable than M200 for constructing the mass function. 

Concerning the mass function of M(< ^^iooo)i it should be reminded that the mass 
function of M200 cannot be analytically derived, and needs to be determined empirically from 
//-body simulations. From this spirit, the mass function of M(< ^1000) can also be obtained 
from iV-body simulations, in a similar manner as the mass function of M2oo- To the best of 
our knowledge, there have not been investigations on the mass function for mass different 
from M200 in literature. Such a new class of mass function, containing less mass scatter 
than the conventional mass function, may be more useful for constraining the cosmological 
parameters. However, whether or not the more accurately measured M{< ^1000) function 
can actually be more useful than the less accurately measured M200 function really depends 
on the detailed form of the M(< ^1000) function that contains the cosmology-parameter- 
sensitive feature, similar to a break at for the M200 function. Therefore our suggestion 
at this point in favor of an alternative mass estimator for a new mass function should be 
regarded as plausible but still preliminary. 

Finally, though the depth of observation has been fixed to i? ~ 25.5 mag in this work for 
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which Ug ~ 30 arcmin"^, for evaluation of mass errors, we may relax this constraint to assess 
other observing depths straightforwardly. At a medium depth R ~ 24.5 mag, suitable for 
wide- field surveys, the background galaxy density Ug ~ 20 arcmin^^ (Fontana et al. 2000). 
One can quickly estimate from Fig. 7 with the Poisson statistics that the galaxy ellipticity 
noise is still much less than the errors at ^^500, ^^looo, and ^1500 introduced by the projection 
effect for clusters of M > 2 x W^Mq and 0.2 < zl < 0.4 and clusters of M > 4 x lO^^M© 
and 0.2 < zl < 0.6. That is, a medium-depth observation can provide as an accurate mass 
measurement as a deep observation for those more massive and lower redshift clusters. Such 
information is crucial for the planning of wide-field surveys. 
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Fig. 1. — An example of the raw k map that produces a negative M200 for the target cluster 
at the center, where the object on the lower-left corner is located at a different redshift. The 
field of view is 1.26'2oo x l-26'2oo defined by the central cluster. 
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Fig. 2. — An example of the raw k map that produces an over-estimated M200; where the 
nearby object in the map is also located at a different redshift. The field of view is defined 
in the same manner as in Fig. 1. 
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Fig. 3. — An example of clean k maps, where contaminants have relatively low signal 
strengths. 
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Fig. 4. — Mass probability distribution functions from noise-free simulations for three dif- 
ferent redshift intervals as indicated in the figure. The error bars represent 68% confidence 
intervals from 20 independent galaxy realizations. 
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Fig. 5. — Mass probability distribution functions of noisy maps for three different redshift 
intervals as indicated in the figure. The error bars represent 68% confidence intervals from 
20 independent galaxy realizations. 
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Fig. 6. — Noise-free [solid) and noisy [dashed) mass probability distribution functions of 
the clean sample for different bins of cluster mass and redshift: 1 < Mi/1O^^M0 < 2, 
2 < Ma/IO^^M© < 4, Ms/W^Mq > 4. Here, the cluster mass refers to Maoo- 
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Fig. 7. — Noise-free [black) and noisy [red) projected- mass radial profiles of the clean 
cluster sample for different mass bins and redshift intervals: 1 < Mi/10^^Mq < 2, 
2 < Ma/lO^^'Mo < 4, Ms/IO^^Mq > 4. Mean values are represented by filled circles (black) 
for the noise-free case, and by connected lines [red) for the noisy case. The error bars indicate 
68% confidence levels around the mean. 
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Tal)le 1: Cluster iiuiiil)er in the raw sample 







0.2<z < 0.4 


0.4 < ^ < 0.6 


0.6 < ^ < 0.8 


1 < Maoo/lO^^AfQ 


< 2 


444 


1144 


1374 


2 < M2oo/1O^^M0 


< 4 


145 


350 


364 


M2oo/10^^Mq > 4 




59 


92 


43 
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Table 2: Cluster number in the clean sample and percentage of interacting clusters among 
those that are removed during the clean 



0.2 < z < 0.4 


0.4 < ;2 < 0.6 


0.6 < ^ < 0.8 


1 ^ i\ J /inl4?\/r oo/( /icO/N 

1 < M200/IU Mq < 2 324 (loyoj 

2 < Maoo/lO^^Mo < 4 91 (20%) 

il72oo/10^M4 > 1 22 (IG'X) 


293 (22%) 

72 (18^) 


38 (25'X) 


Table 3: Dispersion of /i2oo = ^icns(< 6. 


2oo)/Mtruc(< 6'2oo) iu the raw sample 


0.2 < z < 0.4 


0.4 < 2 < 0.6 


0.6 < ^ < 0.8 


1 < Maoo/IO^^M© < 2 (0.69, 0.98) 

2 < M2oo/1O^^M0 < 4 (0.55, 0.72) 
M2oo/1O^^M0 > 4 (0.45, 0.48) 


(0.63,1.13) 
(0.47, 0.76) 
(0.39,0.53) 


(0.74, 1.74) 
(0.65, 1.28) 
(0.36, 0.77) 


For each bin, numbers are given for the noise-free {left) and noisy {right) 


cases. 


Table 4: Dispersion of — -^iens(< ^ 


2oo)/M,rue(< 6*200) in the clean sample 


0.2 < ^ < 0.4 


0.4 < z < 0.6 


0.6 < z < 0.8 


1 < M2qq/W^Mq < 2 (0.60, 0.91) 

2 < M2oo/1O^^M0 < 4 (0.44, 0.63) 
M2oo/1O^^M0 > 4 (0.32, 0.40) 


(0.41,0.68) 
(0.28, 0.42) 


(0.25,0.6) 


For each bin, numbers are given for the noise-free {left) and noisy {right) 


cases. 


Table 5: Mean of /^2oo = -/^iens(< 6*200 )/-^true(< 6*200) in 


the clean sample 


0.2 < z < 0.4 


0.4 < ^ < 0.6 


0.6 < ^ < 0.8 


1 < M2oo/lOi^M0 < 2 (1.02, 1.02) 

2 < Maoo/lOi^M© < 4 (1.00, 0.99) 
Maoo/lO^^Afo > 4 (1.05, 1.04) 


(1.11,1.11) 
(1.08,1.08) 


(1.14,1.12) 


For each bin, numbers are given for the noise-free {left) and noisy {right) 


coses. 


Table 6: Skewness of //200 = Mens(< 6*200 )/-/^true(< 6*200 


) in the clean sample 


0.2 < z < 0.4 


0.4 < z < 0.6 


0.6 < z < 0.8 


1 < Maoo/lOi^M© < 2 (0.6, 0.16) 

2 < M2oo/lO^^Af0 < 4 (0.49, 0.32) 
M2oo/lOi^M0 > 4 (0.60, 0.11) 


(0.56,0.13) 
(0.31,0.02) 


(0.12,0.07) 



For each bin, numbers are given for the noise-free {left) and noisy {right) cases. 
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Table 7: Noise-frcc mean and dispersion for the ratio of the local compensating mass 
AM(6'2oo; l-26'2oo) to the true cluster mass M(< 6200) 





0.2 < ^ < 0.4 


0.4 < ^ < 0.6 


0.6 < ;2 < 0.8 


1 < Maoo/IO^^M© < 2 

2 < M2oo/1O^^M0 < 4 
M2oo/1O^^M0 > 4 


(0.25,0.40) 
(0.21,0.14) 
(0.25, 0.26) 


(0.26,0.16) 
(0.23,0.29) 


(0.18,0.10) 



For each bin, numbers are given for the mean (left) and the dispersion (right). 



Table 8: Mass error dispersions ((t(/X5oo), cr(/xiooo), cr(A*i50o)) in the clean sample for the noisy 
case 





0.2 < 2 < 0.4 


0.4 < z < 0.6 


0.6 < z < 0.8 


1 < Maoo/IO^^MQ < 2 

2 < M2oo/10^^Mq < 4 
M2oo/10^^Mq > 4 


(0.53,0.39,0.32) 
(0.37,0.28,0.22) 
(0.21,0.17,0.13) 


(0.41,0.32,0.28) 
(0.29,0.20,0.18) 


(0.36,0.28,0.25) 



